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Abstract — Since efficient power allocation is a critical require- 
ment for smart grid, we study an important basic setting - 
knapsack problem with selfish users, whereby we design a mecha- 
nism to find a utility-maximizing allocation for a group of users 
with inelastic demands, such that users truthfully reveal their 
private utility information. As a departure from the traditional 
setting, complex-valued entities (e.g. power, voltage, and current) 
are common in A/C electrical grid. There were only few results 
in the literature concerning complex-valued entities for discrete 
optimization, because they are substantially more challenging. In 
this paper, we introduce a non-trivial generalization of knapsack 
problem with a complex-valued constraint on A/C power, which 
casts fundamental insight to discrete optimization for smart 
grid. We provide results of approximability (the existence of 
a (| — e) -approximation algorithm) and inapproximability (the 
absence of FPTAS unless P = NP) for a class of complex-valued 
knapsack problem, considering complex-valued A/C power with 
non-negative real and imaginary parts. Further, we achieve 
truthfulness in this setting. We also apply our results to a setting 
of A/C power with moderate negative imaginary part. 

Index Terms — Complex-valued Knapsack Problem, Discrete 
Optimization for Smart Grid, A/C Electrical Grid, Approxima- 
tion Algorithms, FPTAS, Incentive Compatibility, Truthfulness 

I. Introduction 

The rise of smart grid provides a new breed of resource 
allocation problems. Unlike the traditional settings of resource 
allocation, many entities are expressed by complex values 
in A/C electrical systems, which can conveniently capture 
both phase and magnitude in A/C electrical systems. While 
real-valued entities are often studied in traditional discrete 
optimization, the consideration of complex-valued entities has 
received limited attention, despite its importance to smart grid. 

In this paper, we consider a fundamental problem of utility- 
maximizing allocation of power in A/C electrical grid. Our 
study is based on the well-known knapsack problem El. In 
knapsack problem, there are a number of users; each has an 
inelastic demand of power (expressed by a complex value) and 
a utility value if the demand is satisfied. The goal is to find 
an allocation of power to a subset of users, such that the total 
utility of satisfiable demands is maximized. The constraint is 
that, each inelastic demand is either satisfied or dropped, and 
the total satisfiable demand is at most the capacity of supplies. 

In the terminology of power systems BJ, the real part of 
complex-valued power is known as active power, the imagi- 
nary part is known as reactive power, whereas the magnitude 
is known as apparent power. Electric appliances have different 
active and reactive power requirement^] while the electrical 

'Purely resistive appliances have positive active power and zero reactive 
power. Appliances with capacitive or inductive components have non-zero 
reactive power, depending on the phase lag with the input power. 



grid also has capacity constraints on its active, reactive and 
apparent power supplies. Therefore, the discrete power alloca- 
tion problem in A/C electrical grid is more complicated than 
the traditional knapsack problem, which concerns allocating a 
single real-valued resource. One can formulate discrete power 
allocation problem with separate capacity constraints on active 
and reactive power as two dimensional knapsack problem [4|. 
However, the problem involving capacity constraint on appar- 
ent power appears as a new problem, involving a quadratic 
constraint, rather than linear constraints in traditional knapsack 
problems. In this paper, we provide the first study of complex- 
valued knapsack problem with two important contributions. 

1. (In)approximability: Our first contribution is to provide 
insight on the (in)approximability of complex-valued knapsack 
problem considering active, reactive and apparent power ca- 
pacity constraints. While many discrete optimization problems 
are computationally hard problems for computing the exact 
optimal solutions, approximate solutions that are close to the 
optimal solution may be efficiently computed by a polynomial- 
time algorithm. An effective class of approximation algo- 
rithms is constant approximation algorithms, which produce 
solutions that are within a constant factor from the optimal 
solution in polynomial time of the problem size. Polynomial- 
time approximation scheme (PTAS) can produce a solution 
that is within a factor 1 — e from the optimal solution for 
maximization problems (or 1 + e for minimization problems) 
for any e > within time polynomial in the size of the input. 
The most effective class of approximation algorithms is the 
fully polynomial-time approximation scheme (FPTAS), which 
additionally requires the running time to be polynomial in both 
the input size and 1/e. In this paper, we show the existence 
of a (i — e)-approximation algorithm and the absence of 
FPTAS unless P = NP for a class of complex-valued knapsack 
problems, considering non-negative active and reactive power 
demands. We also apply our results to a setting of A/C power 
with moderate negative reactive power demands. 

2. Truthfulness: Our second contribution is related to 
mechanism design. An algorithm is said to be truthful (or 
incentive-compatible) if none of the users has an incentive 
to misreport his true demand and utility. Namely, they will 
always truthfully reveal their private information. Truthful 
mechanisms are economically desirable in smart grid, because 
it simplifies the design of bidding strategies and facilitates 
users' adaptation. VCG based mechanisms j6j, (7) can realize 
truthfulness, but require exact optimal solutions. However, 
many discrete optimization problems are computationally hard 
(e.g., knapsack problem). Therefore, we often need to adapt an 
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approximation algorithm to achieve truthfulness. Results in the 
literature |1|, [5| give a sufficient condition for truthfulness in 
problem settings with single-minded users: the approximation 
algorithm satisfies a property known as monotonicity . But not 
all approximation algorithms satisfy monotonicity. We show 
that our (| — e)-approximation algorithm for basic complex- 
valued knapsack problem satisfies monotonicity, and hence 
leads to a truthful mechanism. 

II. Problem Formulations and Definitions 



(1) 

(2) 



A. Real-valued One Dimensional Knapsack Problem 
The traditional knapsack problem (1D-KS) is defined as 

(1D-KS) max > XkUk 

subject to ^-^ 

2^ x k d k < C 

where ' 

1) K is a set of users; 

2) ut is the utility of fc-th user if its demand is satisfied; 

3) dk is the positive real-valued demand of fc-th user; 

4) C is a real-valued capacity of total satisfiable demand; 

5) Xk is a decision of allocation determined by an algorithm 
(xk = 1 means that fc-th user's demand will be satisfied, 
otherwise x k = 0). 

This problem can model the discrete optimization of D/C 
electrical systems, where C is the capacity of D/C power. 
1D-KS is NP-complete but has an FPTAS (2), g)- 

B. Two Dimensional Knapsack Problem 

In the setting of A/C electrical systems, a demand for power 
d can be expressed by a complex value d = d n + id 1 . One 
can impose two separate constraints on the capacity of total 
satisfiable demand of active power (i.e. the real part) and 
reactive power (i.e. the imaginary part). This problem can be 
formulated as a two-dimensional knapsack problem (2D-KS): 

(3) 

keK 



(2D-KS) 
subject to 



max 



ax V" x k u k 
{0.1} tTL 



E 

keK 



x k df < C R and ^ x k d\ < C 1 



(4) 
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2D-KS can be generalized to m-dimensional knapsack 
problem (mD-KS). toD-KS is NP-complete since 1D-KS is 
a special case. mD-KS only has a PTAS, and there is no 
FPTAS for toD-KS, when m > 2 (unless P = NP) g). 

C. Complex-valued Knapsack Problem 

Our study concerns the capacity on apparent power (i.e. the 
magnitude of the total satisfiable demand). We formulate a 
complex-valued knapsack problem (C-KS) as follows: 

(C-KS) max } x k u k (5) 

subject to | I 

x k d k \ < C (6) 



keK 



where 

1) dk = rffe + is the complex-valued demand of power 
for fc-th user; 

2) C is a real-valued capacity of total satisfiable demand 
in apparent power. 



Evidently, C-KS is also NP-complete, since 1D-KS is a 
special case when we set all d\ = 0. 

We also define a general setting of knapsack problem 
for A/C electrical systems with a combination of capacity 
constraints on apparent power, active power and reactive power 
together. Namely, we formulate a general complex-valued 
knapsack problem (GC-KS) as follows: 



(GC-KS) 
subject to 

^ x k d k 



max 

x fc e{o,i} 



keK 



XkUk 



(7) 



keK 



< C and J2 x k d f < c * and J2 Xkd k - ^ 

keK 



keK 



(8) 

One of our contributions is to provide an approximation 
algorithm for C-KS and GC-KS. 

D. Approximation Algorithm 

Given an allocation {x k }, we equivalently represent it by 
the satisfied subset of users S C K, where S = {k £ K : 
x k = !}■ For a subset S, we denote u(S) = J2kes Uk - 

We denote S* to be an optimal solution of C-KS. An algo- 
rithm A for C-KS is said to be a p-approximation for some 
constant p(< 1), if S is the output of A and u(S) > p-u(S*) 
on each input. 

For the convenience of analysis, in this paper we assume 
that each demand d k lies in the first quadrant, such that 



df > and d\ > 



(9) 



We then provide a (i — e)-approximation algorithm to solve 
C-KS under this assumption. 

Note that our approximation algorithm is also applicable to 
the setting of negative reactive power, such that for all fc 



d% > and d% > \d\\ 



(10) 



The basic idea is to rotate the first quadrant by an angle of f 
clockwise. Normally, positive active power demand implies 
a power consumer (who requires positive power from the 
power grid), while negative active power demand implies a 
power supplier. Hence, it is often that d R > 0. For reactive 
power, inductors have positive reactive power, while capacitors 
have negative reactive power. In most electrical power systems 
regulations, it is often required that the absolute reactive power 
is not too large, as compared to active power [3|. Thus, the 



assumption in Eqn. (IOi is justified in common settings. 



E. Truthfulness 

Suppose the set {(dk,Uk)} are users' private information. 
Without the precise knowledge of the private input, an al- 
gorithm cannot find a true utility-maximizing allocation. We 
seek to design a truthful mechanism in a way that incentivizes 
all users to report their true private information. A general 
approach is to enforce payment on each user that depends 
on reported inputs, such that revealing the true information 
maximizes his net utilit^j^] (utility minus the payment). In this 
way, no user has an incentive to misreport his private [dk, Ufe). 

2 Here, to conform to the traditional usage of utility in knapsack problems, 
we change several important names in mechanism design: our net utility is 
the utility in mechanism design, and our utility is called valuation. 
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There is rich literature on designing truthful (or incentive 
compatible) mechanisms for problem settings well studied in 
traditional algorithmic design. The setting in C-KS is a special 
case of single-minded users. (See Q], (5J) 

For each user k G K, we let tk = (dfe,itfc). Let t = {tk} 
and t- k = {(d 3 ,Uj) : j G #\{fc}}. 

We define a mechanism as Ai — (A,p), consisting of an 
algorithm A(t) that computes an allocation {xk} and payment 
Pk(t) for each user k. We write -4fc(i) = x k - 

Since is private information, user k may not truthfully 
reveal tk- For any input of other users t-k, the utility of user 
k depends on his report t' k = (d' k , u^): if Ak{t) = 0, he is not 
satisfied and utility equals 0; if d' k < d^ or d' k < d\, i.e., 
the power, if allocated to him, would not be enough for his 
demand, his utility is also 0; otherwise, he gets Uk- 

Clearly, truth-telling results in net utility Uk(t) — Ak(t) ■ 
Uk —Pk(t)> if user k misreports t' k , his net utility Uk(t' k , t-k) 
becomes the utility defined above minus £>fc(ij.,t-fc). 

A truthful mechanism maximizes user k's net utility when 
he reports true tk for every t-k, tk and t' k : 



U k (t k ,t-k) > U k {t'k,t- k ) 



(11) 



An algorithm A is said to be truthful, if there exists 
payments p such that mechanism M. = (A,p) is truthful. 

A sufficient condition to ensure truthfulness for single- 
minded users is monotonicity Q], j5). In our setting, 
monotonicity is specified as follows. 

Definition 1: An algorithm A is monotone if Ak(t) = 1 
implies A k (t' k ,t- k ) 
d'\ < d\. 



1 for any u' k > u k , d' k < d k 



and 



Intuitively, in a monotone algorithm, if user k is satisfied with 
demand dk and utility Uk, then he should be also satisfied 
when he has smaller demand and larger utility. The following 
theorem states the sufficiency of monotonicity QJ, Q. 

Theorem 1: (T], (3} Let A be a monotone algorithm for 
single-minded users. Then there exists payment function p 
such that M — (A,p) is truthful. 

The basic idea is that, given a monotone algorithm A, if 
we fix dk and t-k for user k, then A defines a critical value 
9k(dk,t-k), such that when Uk is above the critical value, k 
is satisfied; and when Uk is below the critical value, k is not 
satisfied. Then we can define a payment function pit), such 
that each satisfied user pays the critical value: 

' k (d k ,t-k) if Ak{t) = 1 
otherwise 



Pk(t) 



(12) 



By Theorem [T] if we are able to design a monotone 
algorithm, we can transform it into a truthful mechanism. 
Moreover, the critical value for a given input can be computed 
in polynomial time by a binary search on interval [0, Uk] for 
each user k during which we repeatedly test if k is satisfied by 
running algorithm A. Therefore, a monotone polynomial time 
algorithm A implies a polynomial time truthful mechanism. 



While monotonicity usually holds for exact optimization, 
many well-known approximation algorithms violate mono- 
tonicity. One example is the FPTAS for 1D-KS B), which 
is based on scaling utilities and applying the optimal pseu- 
dopolynomial dynamic programming algorithm, but is shown 
to be not monotone. Recently, (TJ adapts it to a sophisticated 
scheme, which is a monotone FPTAS for 1D-KS. 

III. Monotone {\ - e) -Approximation Algorithm 

We present a polynomial -time monotone ^-approximation 
algorithm for C-KS, which relies on a polynomial-time mono- 
tone approximation algorithm for 1D-KS as a subroutine. 

A. Basic Idea 

We first present the intuition of our algorithm. Recall that 
we assume d^ > and d\ > 0. Since we may preprocess 
the demands and eliminate every demand whose magnitude 
exceeds capacity C itself, without loss of generality, we 
assume all \d k \ < C. 

Graphically, each demand is a vector d = d R + id 1 in the 
first quadrant. A feasible solution of our problem is a subset 
of users whose sum of demands lies in region V, the 1/4 disk 
of radius C in the first quadrant. As shown in Fig. [T] T> is 
divided by chord PQ into triangle T)\ and circular segment 
T>2- The ? line intersects chord PQ at point R. 




Fig. 1. An illustration of notations. 

Our algorithm is based on the following two observations: 
1) If we project all demands onto the ? line, i.e., 



(df + dl 



(13) 



we make all demands 1 -dimensional. Now a subset of 
demands has sum 

E dk < C/V2 (i.e., the sum vector 
does not go beyond point R on the j line) if and only 
if its original sum vector E dk lies inside the triangle 
T>\. This is because that, the sum of projections, E dk, 
is the projection of E dk on the ? line. Therefore, the 
subproblem on feasible region T)\ can be solved by an 
algorithm for 1D-KS with demands changed to dk and 
capacity to C/V2. 
2) The subproblem on feasible region T>\ is almost the 
whole story: First, an optimal solution in T> can contain 
at most one demand in T>2', second, if an optimal solution 
consists of more than one demand, its sum can be 
broken into either two separate subsums lying in T>x, 
or, the sum of a vector in T>2 and a subsum in T>\. 
Our algorithm essentially approximates the optimum 
between an optimal solution for the subproblem on 
feasible region T>\ and an optimal solution on input 
demands lying in T> 2 - This gives an approximation ratio 
of 1/2- e. 
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B. Approximation Algorithm 

We let Alg a [(d fc : k £ K),C] be our algorithm for C-KS, 
where (dk ■ k £ K) is the input of the demands and C is the 
capacity. Also, we let Alg ld [(<ifc : k £ K), C] be a polynomial- 
time approximation algorithm for 1D-KS (for example, from 
Q). We describe our algorithm as follows: 

Algorithm 1 Alg a [(4 :k£K),C] 



for k £ K do 

Set dk = min{ 
end for 

Set 5 = Alg ld [(4 
Output S 



4+4 ci 



In Alg a , we first project all demands onto the j line. If the 
projection of a demand is larger than C/y/2, we cut it off to 
C/y/2. Then we use an approximation algorithm Alg ld for 
1D-KS to compute an allocation considering the projected 
demands and capacity C/y/2. 

Theorem 2: If Alg ld is a p-approximation algorithm for 
1D-KS, then Alg a is a ^-approximation algorithm for C-KS. 

Corollary 3: Since 1D-KS has an FPTAS (TJ, there is 
a (| — e) -approximation algorithm for C-KS that runs in 
polynomial-time in the size of input and 1 /e, for any e > 0. 

C. Proof of Theorem [2] 

Let S be the output of Alg a and S* be an optimal solution 
to C-KS, for which the feasible region is T>. We partition K 
into two disjoint sets K\ and K 2 , such that K\ = {k £ K : 
d k £ £>i} and K 2 = {k £ K : d k £ V 2 }. Note that the 
projection of any demand in K\ onto the f line is at most 
C/y/2, whereas that in K 2 is larger than C/y/2. 

Let Si be the output of Alg a , when the input is K\. Let 
S2 be the output of Alg a , when the input is K 2 . Let S* 
and S 2 be their corresponding optimal solutions. SJ is an 
optimal solution to the 1D-KS on projected dem ands w ithin 
capacity C/y/2, hence by our observations in Sec. 
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optimal solution to the C-KS on feasible region T>i. On the 
other hand, since each demand in K 2 is changed to a vector 
exactly equivalent to the capacity limit of the 1D-KS, only 
one of them can be satisfied. Hence S 2 chooses the one with 
maximum utility S 2 = {arg max^g^ u k }. 

Since any demand in K 2 will not combine with any in K\ 
to form new feasible solutions to the 1D-KS, Alg a outputs 
either a solution whose sum vector lies in T)\ or a singleton 
set of a demand in K 2 , which is evidently a feasible solution 
to the C-KS. 

Optimally the 1D-KS outputs arg max{u(Si ), u(S 2 )}. 
Since Alg ld is a p-approximation algorithm to the 1D-KS, 
we have u(S) > p ■ max{u(Si ), u(S 2 )}. 

Next, we analyze the approximation ratio of Alg a in three 
cases. Here for a subset S C K, we denote 

fees fees fees 
Case (1): (p-approximation) We consider an optimal solution 
S*, such that the sum of demands d(S*) £ X>i. 



Proof: This is an easy case u(S*) = u(S*). We have 
u(S) >p-max{u(St),u(S%)} >p-u(Sf) = p-u(S*). ■ 
Case (2): (^-approximation) We consider an optimal solu- 
tion S*, such that d(S*) £ T> 2 , and there exists a user j £ S* 
whose demand dj £ T> 2 . 

Proof: Let z = 2~2keS*/{j} dk- Thus, d(S*) = dj + z. 
That is, the sum of demands of S* can be written as the sum 
of demand dj and a demand subset sum z. (It is possible that 
S* only consists of a single user j.) Note that dj £ T> 2 and 
z € T>\. Otherwise, the projection of d(S*) = dj + z on the 
I line would exceed 2 • C/y/2 > C. 

Moreover, we have u(S*/{j}) < u(S*), because SJ is an 
optimal solution for feasible region T)\. On the other hand, 
Uj < u(S|) by S 2 — {argmaxfe e i<- 2 Uk}- We obtain: 

u(S*) = Uj + u(S*/{j}) < u(S* 2 ) + u(St) (15) 

Combining with max{M(Si ), u(S 2 )} > \ (u(S*) + u(S 2 )) 
and u(S) > p ■ max{u(S?), u(S%)}, we get u(S) > fu(S*). 

■ 

Case (3): (^-approximation) We consider an optimal solu- 
tion S*, such that d(S*) £ T> 2 , and dk £ T>\ for every user 
k£S*. 

Proof: First, we define d(S) = 2~2kes^ k - ^ projection 
on the j line, it is equivalent to saying that d(S*) > C/y/2, 
and d k < C/y/2 for every user k £ b . 

The proof of Case (3) is immediate from Lemma [T] as 
follows: By Lemma[l] we have d(T) < C/y/2 and d(S*/T) < 
C/y/2 for some subset T C S*. That is, d(T) £ V x and 
d{S*/T) £ V x . 

Thus, u(T) < u(Sf) and u(S* /T) < u(S$). Moreover, 
u(S*) = u(T)+u(S*/T). Hence, u(S*) < 2it(S|). 

Finally, the total utility of the output solution of Alg a , 

u(S) > p-max{u(S 1 *),w(S 2 *)} >p-u(&t) > ^u(S*) (16) 

■ 

Combining Cases (l)-(3): min{p, p/(l + p),p/2} = p/2, 
we complete the proof of the approximation ratio of Alg a as 
p/2. 

Lemma 1: For a set of n positive real- valued numbers 
ai,...,a n satisfying 2~27=i ai — a>i < C for all i and 
C > C/y/2, there exists a subset T C {1, ...,n} such that 

^a, <C" and ^ a, < C (17) 

ieT ie{l,-,n}\T 

Proof: Let j be the smallest index such that the partial 
sum exceeds C , i.e., X^=i a i — C an( ^ 2~2i=i a i > C' . 
Clearly j > 2 since all dj < C. 



and y = E"=,-+i a i 



Let x = Y?i=i a i> z z 

Note that Y^i=i ai — x + V + z - We already have 

x < C, z<C, x + y + z>C' and x + z > C' (18) 

The lemma holds if y + z < C, because we can set T = 

If y + z > C, then we obtain: 

x + y = 2{x + y + z) - (x + z) - (y + z) (19) 

< 2C- 2C < (2 - y/2)C <^< c ' ( 20 ) 
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because x + y + z < C. Hence, we can set T = {1, j — 
l,j + l,...,n}. ■ 

D. Monotonicity 

We also show that our algorithm Alg a satisfies monotonicity. 

Theorem 4: There is a truthful (1/2 — e)-approximation 
algorithm for C-KS that runs in polynomial time in the size 
of the input and 1/e, for any e > 0. 

Proof: Note that 1D-KS has a monotone FPTAS from 
JlJ. We use it as Alg ld . Theorem [4] follows immediately from 
Theorem[T]and Corollary[3] if we prove that Alg a is monotone. 

We then need to show that, if user k is satisfied by Alg a 
with demand dk and utility Uk, k is also satisfied with demand 
d' k and utility u' k , where u' k > Uk and d'k < d^ and d'\ < d\, 
while all inputs of other users do not change. 

User k is satisfied by Alg a on d' k and u' k if and only if it is 

satisfied by Alg ld on d' k and u' k . Since d k = min{^-±^, 

d'^ < dfc and d'\ < d\ implies that d' k < d k . Then from the 
monotonicity of Alg ld , k is satisfied by Alg ld , and hence by 
Alg a . ■ 

IV. INAPPROXIM ABILITY 

In this section, we complete the study of C-KS by providing 
an inapproximability result. We show that C-KS can not have 
FPTAS, unless P = NP 

We remark that it is known there is no FPTAS for 2D-KS 
(see [4]), which does not have direct implications for C-KS. 
Neither can the reduction in its proof be directly adapted. Our 
proof is a clever modification of the reduction for 2D-KS. 

We reduce the EQUIPARTITION problem to C-KS: 

Definition 2: (EQUIPARTITION Problem): Given a set of 
positive integers {u>k ■ k G K}, with \K\ = n even, we 
determine if there is a subset of items S C K such that 



i 5 i = 2 and J2 Wk 

kes 



k<£S 

It is well-known that EQUIPARTITION is NP-complete. 



Theorem 5: There is no FPTAS for C-KS, unless P = NP 
Proof: We define a decision version of C-KS with a 
cardinality objective: given {wk ■ k G K}, a capacity bound C 
and a cardinality bound M, we determine if there is a subset 
of items S such that 



\S\ > M, and 



fees 



d k < C 



Now we map every instance of EQUIPARTITION to an 
instance of the C-KS decision problem that always yields the 
same answer. 

Given {u> fc : k £ K} from EQUIPARTITION, define 



M = n/2, 



Wk, 



P(w D 



C = 



v(t) +/? I" 



w 
Y 



where W = YJl=i w k, w n 
in our reduction, d\ > 0. 



= maxjwfe : k G K}. Note that 



As shown in Fig. [2] the feasible region T> for C-KS is the 
2 disk of radius C in the first quadrant. Since for any subset 

s c k, J2 kes 4 = /W-^nuBt-Efces the cardinalit y 

constraint \S\ > f imposes all solutions to have its sum vector 
in the halfplane H : d 1 > /3(^f^ - d R ). The dividing line 

of H goes through point P : (^-,/3(^f^ - -f )). Our main 
idea is to set f3 > such that the dividing line of H coincides 
with the tangent line at P. Thus we make the intersection of H 
and T> exactly P, which implies \S\ = ^ and 2~2kes w j = "T 
for any solution S to our reduced C-KS decision problem 
instance. 




Fig. 2. Reduction of inapproximability. 



On the other hand, it is clear that each subset S satisfying 
conditions of EQUIPARTITION also satisfies conditions of the 
reduced C-KS decision problem. Therefore, the solution of the 
reduced C-KS decision problem is equivalent to the solution 

of EQUIPARTITION. 

To determine a proper f3, since the dividing line of halfplane 
H goes through P, it coincides with the tangent line at P if 
and only if they have the same slope, i.e., 

w 

2 _ a 



PC- 



x) max _ W 
2 2 



Solving the above equation, we obtain 



P 



W 



W 



which is > unless all weights are equal. In this case, we set 
P = 0. It is trivial to check that both instances have answer 
"Yes". 

So far we have shown the NP-completeness of the complex 
valued knapsack decision problem. So its maximization ver- 
sion, where \S\ > M is replaced by max is NP-hard. This 
is equivalent to the original complex-valued problem with all 
Uk = 1. 

Finally, we use the standard technique to prove the inap- 
proximability by FPTAS (for the above maximization version). 
Suppose that there exists an FPTAS for any e > in time 
polynominal in n and 1/e. Then we choose e = ^q-p Let the 
sum of values of the optimal solution be z* > and that of 
the approximation solution produced by FPTAS be z A . We 
obtain: 



> (1 — e)z* > z* — z* jri > z* 



1 



because z* < n. Moreover, since z* is an integer, this implies 
that the FPTAS can solve the problem exactly in polynominal 
time, contradicting the NP-hardness of the problem. Therefore, 
there is no FPTAS to C-KS. ■ 
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V. Generalized Complex- valued Knapsack 

We are also able to solve the generalized problem GC-KS 
by changing our approximation algorithm Alg a in Section III 
Now, instead of an approximation algorithm Alg ld for 1D- 
KS as a subroutine, we rely on an approximation algorithm for 
3D-KS (three dimensional knapsack problem) as a subroutine. 
3D-KS is defined as: 

(3D-KS) max V x k u k (21) 

subject to keK 
£ x k d\ < C 1 and ^ x k d\ < C 2 and ^ x k d\ < C 3 

keK keK keK 

(22) 

where d k > for j — 1, 2, 3. 
A. Basic Idea 

Now a feasible solution of our problem is a subset of users 
whose sum of demands lies in the intersection of halfplanes 
d R < C R , d l < C\ and the 1/4 disk of radius C in the 
first quadrant. In the most general case (C R , C 1 < C), both 
halfplanes cut the circle, which also cut the original regions 
T>i and T>2 in Section III Fig. [3] shows the new T>i (polygon 
PSTQO) and T>2- Clearly, the feasible region T> is the union 
of T>i and T> 2 . 




Fig. 3. An illustration of notations. 

Recall that OR is perpendicular to ST and the length of OR 
is C/^/2. If we denote the projection of a demand d k onto line 
OR by d k , the region T>i corresponds to 3-dimensional linear 
constraint J2keK x kd~k < C/y/2, J2keK x k d k < C R and 
TlkeK x kd\ < C l . Thus the subproblem on feasible region 
T>i can be solved by a 3-dimensional knapsack algorithm. 

On the other hand, the solutions in polygon T)\ is almost 



the whole story by the same reason as in Section III Here 
our algorithm directly takes the maximum between an optimal 
solution for the subproblem on feasible region T>\ and an 
optimal solution on input demands lying in 2?^] This gives 
an approximation ratio of 1/2 — e. 

The degenerate cases (C R > C or C 1 > C or both) can be 
treated easily by setting T, Q to be the intersection point of the 
circle and the d R -axis, or setting P, S to be the intersection 
point of the circle and the d^axis, or both. 

B. Approximation Algorithm 

Let Alg c [(d fe : k <E K),C,C R ,C I ] be our approxima- 
tion algorithm for GC-KS. Let Alg 3d [((4, d 2 k , d%) : k £ 
K),C X ,C 2 , C 3 ] be an approximation algorithm for 3D-KS 
(e.g., from [2] or [4]). We describe our approximation algo- 
rithm to GC-KS by Alg c . 

3 We can also modify Algorithm 1 of Section |lll| in this way, which is 
more intuitive and improves the ratio of case (2) to p/(l + p), but has no 
monotonicity guarantee. Since we are not aware of any monotone PTAS for 
3D-KS, we take this alternative here. 



Algorithm 2 A\g c [(d k 
l: for k e K do 

Set d k = 



k G K),C,C R ,C 1 } 



end for 

Set Si = Alg 3d [((J fe ,d R 4) : k e K), ^C R ,C l ] 
Set = {a,Tgmax keK . dkeT)2 {u k }} 
Set S = argmaxs l! s 2 {u(S , i),u(S , 2)} 
Output S 



Theorem 6: If Alg 3d is a p-approximation algorithm for 
3D-KS, Alg c is a ^-approximation algorithm for GC-KS. 

Corollary 7: Since 3D-KS has a PTAS g), there is a 
— e) -approximation algorithm for GC-KS that runs in 
polynomial-time in the size of input, for any e > 0. 

The proof of Theorem [6] is similar to that of C-KS. We 
omit the proof here because of the space limit. Note that 
here u(S2) = v-iS^) since 5*2 only consists of a single 
demand in T>2 that has maximum utility. And u(S) = 
max{u(5'i), u(5 < 2)} > max{p • u(Sl), u(S*2 )}■ 
VI. Conclusion 

Balancing the supplies and demands in smart grid becomes 
increasingly challenging, which highlights the importance of 
our fundamental study of discrete resource allocation in A/C 
electrical systems. We investigate a basic problem underpin- 
ning smart grid operations, complex-valued knapsack problem, 
a departure from the traditional real- valued resource allocation 
settings. Our results of (in)approximability and truthfulness 
provide the basic insights on efficient and economically vi- 
able mechanisms for discrete optimization in A/C electrical 
systems. We envision a wide range of applications that build 
on our results, such as (1) scheduling of EV charging with 
inelastic commuting distance constraints, (2) scheduling of 
heavy industry productions that require large minimum power 
supplies. In the future, we will work on the extension of our 
results to a networked setting, which is a generalized integer 
multi-commodity problem with complex-valued supplies and 
demands. We also seek to find a PTAS for our problem, closing 
the gap between constant approximation and FPTAS. 
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